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Searches for SUSY Particles at HERA 

Y. Sirois a * 

a LPNHE, Ecole Polytechnique, IN2P3-CNRS, 
91128 Palaiseau, France 

Searches at the ep collider HERA for supersymmetric partners of the Standard Model fermions are presented. 
Assuming R-parity conservation, selectrons and squarks of the Minimal Supersymmetric Standard Model are 
excluded for masses up to 65 GeV in a new region of the standard parameter space. Admitting in addition a 
R-parity violating Yukawa coupling A', squarks are sought through direct e-q fusion in a yet unexplored mass- 
coupling domain, taking into account all possible squark decays. Squark masses up to 240 GeV are excluded for 
A' > \/4vraem, depending only weakly on the free model parameters. 



1. INTRODUCTION 

Supersymmetry (SUSY) is a likely ingredient 
of a true fundamental theory beyond the Stan- 
dard Model (SM). Until the SM gets falsified by 
experimental data, the Minimal Supersymmetric 
extension of the Standard Model (MSSM) offers 
a natural framework to guide the search for new 
scalar superpartners of SM fermions. 

The analyses described in this paper assume 
the minimal field representation of the MSSM. 
This contains in particular selectrons ez, and en, 
and squarks qx and Qr, as partners of left- and 
right-handed electrons and quarks. The masses 
of such bosonic sparticles are treated here as free 
parameters. Gauginos and higgsinos, respectively 
partners of the Z°, 7, W /± gauge bosons and 
of the two required Higgs doublets, mix in the 
neutral and charged sectors to form four neu- 
tral X? 2 3 4 (neutralinos) and two charged X12 
(charginos) mass eigenstates. The mass mixing 
matrices for the Xm an d Xn depend on the soft 
supersymmetry breaking parameters Mi and M2 
associated to SU2 and U\ gauginos, on the ra- 
tio tan/3 = V2/V1 of the two Higgs field vac- 
uum expectation values, and finally on the mix- 
ing parameter fi associated to Higgs superfields. 
In the strict MSSM context, a discrete sym- 
metry ensures the multiplicative conservation of 
the -R-parity quantum number defined as R p = 
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(— 1) 3B+L+2S , where B denotes the baryon num- 
ber, L the lepton number and S the spin such 
that Rp = 1 for the SM particles and —1 for 
their superpartners. As a consequence, sparti- 
cles must be produced by pairs and the lightest 
SUSY particle (LSP) is stable. Throughout this 
paper, it is assumed that the LSP is the lightest 
neutralino Xi- The number of free MSSM param- 
eters is further reduced by assuming that Mi and 
M2 are related via the weak mixing angle 9w as 
suggested by Grand Unified Theories (GUT's), 
Mi = 5/3Af 2 tan 2 9 W . No other GUT relations 
are used. In particular, the mass of the gluinos 
g, partners of SM gluons, is left here as a free 
parameter. We assume that M§ > Mq such that 
q —t qg decays are kinematically forbidden. 

The general SUSY superpotential allows, be- 
yond the MSSM, for gauge invariant terms with 
R p violating (fi p ) Yukawa couplings between the 
quarks and leptons and their squark (g) or slep- 
ton (I) partners. Of particular interest for HERA 
are the terms X' i; j k LiQjDk which allow for lep- 
ton number violating processes. Provided that 
baryon number is effectively conserved at low en- 
ergy, sizeable A' couplings are consistent with 
GUT's, Supergravity and Superstring theories 0. 
They moreover do not unavoidably suffer from 
cosmological constraints ^] and are even required 
for baryon asymmetry genesis in some cosmo- 
logical models with first order electroweak phase 
transition By convention the ijk indices 

correspond to the generations of the superfields 
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Li, Qj and Dk containing respectively the left- 
handed lepton doublet, quark doublet and the 
right-handed quark singlet. Expanded in terms 
of matter fields, the corresponding interaction La- 
grangian is B : 
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LiQjD k 



= A' 



ijk 
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-v L d J L d k R 



h.c. 



where the superscripts c denote the charge con- 
jugate spinors and the * the complex conjugate 
of scalar fields. Hence the couplings Ay fc allow 
for resonant production of squarks through e-q 
fusion. 

Slepton-squark searches at HERA in the 
MSSM framework are described in section ^. 
Squark searches in the $ p extension of the MSSM 
are discussed in section ||. Conclusions are pre- 
sented in section [| 

2. SEARCH FOR SLEPTON-SQUARK 
PRODUCTION IN THE MSSM 

The dominant MSSM process at HERA is the 
pair production ep — ► eqX via i-channel \m ex ~ 
change as shown in Fig. [I]. The production cross 
section || depends on first approximation only on 
the Xi mass and couplings, and on the hard pro- 
cess energy threshold (Mg + Mg). The sfermion 
decays involve mainly the \i which lead to strik- 
ing event signatures looking like Deep Inelastic 
Scattering Neutral Current (DIS NC) except for 




Figure 1 . Lowest order Feynman diagram for e + q 
production in ep collisions with typical subsequent 
sfermion decays involving the LSP. 



large associated missing energy and momentum. 

The HI analysis || based on an integrated lu- 
minosity of L = 6.4 pb -1 , requires an isolated 
e within 10° < 8 e < 135° possessing an en- 
ergy E e > 10 GeV and a transverse momentum 
P± > 8 GeV. Furthermore, it is asked that the 
quantity — P z which peaks at twice the in- 
cident electron energy (i.e. 55 GeV) for DIS NC 
events be found below 40 GeV. Finally, the abso- 
lute value of both the parallel and perpendicular 
(relative to the measured electron) component of 
the total transverse "momentum" P± should be 
large (> 3 7 GeV). 

No event candidate remains in the data while 
0.6±0.2 radiative charge current (CC) events that 



> 
a) 



60 
?F 50 
40 
30 
20 
10 




tanp = 1 .41 




H1 


100 > 

0) 

C3 








80 5 ~ 








60 




|r < -50 GeV A 




40 












IS 
IT 

r / 

Ill 


? - 

a. 


20 



30 



40 



50 



60 



70 



(M i + M-)/2 (GeV) 



> 



100 




-300 



-200 



-100 
H (GeV) 



Figure 2. Excluded mass domains at 95% CL in the 
(top) M x o versus (Mg + Mq)/2 plane and (bottom) in 
the M2 versus fi plane for tan/3 = 1.41. 
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mimic the desired event topology are expected. 
This background dominates over expectations 
from DIS NC and j-g fusion (photoproduction) 
processes. The signal detection efficiency e(V) 
with V 2 = (M? - M 2 „ ) (M| - M 2 ) /4M s M q is es- 

Xi y Xi 

timated to reach 56% for V > 25 GeV and drops 
to < 20% for V < 10 GeV. 

Exclusion limits are derived at 95% Confidence 
Level (CL) and presented in Fig. |[ The excluded 
mass range (Fig. | top) extends to (Mg+M§)/2 = 
65 GeV and to M v o = 40 GeV. The search is sen- 

Xi 

sitive down to mass differences of M x o — Me } q ~ 
10 GeV. In the MSSM parameter space, the ex- 
cluded domain (Fig. || bottom) covers an area at 
/i«0 where the Xi 1S dominated by its 7 com- 
ponent such that the couplings of the e and q 
are of electromagnetic strength and the cross sec- 
tion is maximal. The domain excluded by HI 
results extends to larger values of M% than recent 
MSSM searches at LEP results @] and might still 
be uncovered by searches at the TEVATRON for 
Mg > Mg and M x o ~ M g @. 

3. SEARCH FOR # p -SUSY SQUARKS 

Squarks can be directly produced in the s- 
channel at HERA through e-q fusion as shown 
in Fig. ||. By gauge symmetry, only partners of 
left-handed it-like quarks (i.e. ul, cl or t^) or 
right-handed d-like quarks (i.e. cLr, sr or 6r) 
can be resonantly produced through nine possi- 
ble X[j k couplings. The production cross-section 
scales with A' 2 and with the probability q'{x) to 
find in the proton the relevant quark with mo- 
mentum fraction x ~ M~/s ep for ep collisions at 
^fs^ — 300 GeV. The couplings A' llfc are best 
probed with an e~ beam via processes involving 
the u valence quark (e.g. e~ + u — > Sr) while 
Ay j (j 7^ 1) are best probed with an e + beam 
via processes involving the d valence quark (e.g. 
e + + d — » cl). The other A' lj7c couplings (j, k ^ 1) 
involve heavier quark flavours that have to be ex- 
tracted from the sea. 

The dominant squark decays result from a com- 
petition between the $ p Yukawa couplings and 
the MSSM gauge couplings. "Right"- and "left"- 
handed squarks have different allowed or dom- 
inant decay modes. In "^ p -decay modes", the 



ML-like squarks decay only via ill — ► e + d while 
the d^-like squarks have an additional decay with 
missing transverse momentum due to an escap- 
ing neutrino oIr — > e + u or vd. The partial de- 
cay widths for these $ p -decay modes scale with 
A' 2 x Mq. In "gauge decay" modes, the UL-like de- 
cays involve either the \m or the %+, % — ► uXm 
or dXn- ln contrast the gauge symmetry forbids 
an electroweak coupling of the dp-like squarks to 
the xt- Hence the cjr — ► q'xt ^ s strongly sup- 
pressed for first and second generation squarks 
as it can only proceed through the H + compo- 
nent of the x + with a coupling proportional to 
the q' mass. The partial decay widths for the 
gauge decays scale with MqX (1— M 2 /M?) 2 and, 
in detail, depend on MSSM parameters via the 
mass- mixing in the gaugino-higgsino sector 
The gauge decay modes are expected to domi- 
nate through most of the parameter space for the 
A' — Mq domain accessible at HERA. The fy p - 
decay modes will only dominate near the kine- 
matical limit (e.g. at Mq > 250 GeV) where a 
large A' is required at production. 

The ul gauge decays dominantly involve a xt 
as soon as kinematically allowed, i.e. everywhere 
except in a domain at large M2 and /i«0 where 
the Xm i s heavy and has a large 7 component. 
This can be seen in Fig. [|. The lightest state 
Xi is found to be concerned in a vast portion 
of the parameter space. The dR gauge decay in- 
volves dominantly the x? (> 90% of gauge de- 
cays) when this LSP is dominantly 7-like, i.e. at 
(J-/M 2 < —5/9. Else the LSP participates in 20 to 
80% of the gauge decays in the region where its 
Z component dominates. Its contribution is neg- 
ligible only in the /i < region when dominated 
by its H component. 

The x' s couple to fermion-boson pairs. The 
X + can decay into a x° an d two matter fermions 
(quark or lepton pairs) via the exchange of a real 
or virtual W + boson, e.g. xt ~ * y?[W + \W + — > 
qq' or l + v. It can also decay into three matter 
fermions (a lepton and two quarks) when involv- 
ing the JjL v couplings Ay & through the exchange 

of a virtual sfermion, e.g. xt ~ * du^UL — > e + d 

01 xt ~ * u d~L ; d~L — > vd. The $p decays xt ~ ¥ 
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Figure 3. Lowest order s-channel diagrams for ul 
squark production at HERA followed by (top) a Ijt p 
decay or (bottom) a gauge decay involving a gaugino- 
higgsino x- Typical decays of the emerging Xi or xt 
are shown in the dotted boxes. 
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Fi gure 4. (top) Xm or Xn eigenstates involved in 
the dominant decay modes of the ul in regions of the 
(M2,/j) plane for tan/3 = 1; (bottom) physical masses 
of the Xm or Xn versus fu, for M2 = 60 GeV. 



e + djdk or v e Ujdk dominate for Ayj. couplings val- 
ues of 0(-J 'An The Xi is generally unsta- 
ble and will also decay via the IjL p coupling X[j k 
through virtual sfermion exchange, X? ~ * vdjdk, 
e + iijdk or e~Ujdk- The X? being a Majorana par- 
ticle, there is an equal probability for decays into 
the like (with respect to the incident e beam) and 
unlike sign leptons. Thus, this leads to striking 
and background free lepton number violating sig- 
nals. The Xi - * e^qq decays dominate ( 63% to 
88% of the x? branching ratio ) where the x? is 
7-like. The X? - * vqq decays generally dominate 
where the x? is Z-\ike. A //-like x? w iU most 
probably be long-lived and escape detection. 
In summary, the event topologies can be 



broadly classified in three families. (A) squark IjL p 
two-body decays (e.g. q' — ► vq' or q' — ► eq); (B) 
squark gauge decays into three charged fermions 
final states (e.g. q' — > q'XilXi ~* e ± q"q"') or for 
{ti-like squarks, (e.g. <f -> qxf;Xi ~> e+ 9?) ( c ) 
squark gauge decays with escaping u's (missing 
transverse momentum p±), jet(s) and/or charged 
lepton(s). These were studied in a recent HI 
analysis based on L = 2.8 pb _1 . The event 
topologies of type (A) are indistinguishable on 
event-by-event from the DIS NC or CC back- 
grounds. These are statistically suppressed by 
requiring either an isolated e with > 7 GeV in 
the angular range 10° < 8 e < 145° (e + q + X final 
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states) or and overall^ > 25 GeV (v+q+X final 
states) . In the former case, the background is fur- 
ther suppressed by requiring that the e be found 
at high y e where y e = 1/2(1 + cos 0*) and 9* e is 
the electron angle in the e-q CM frame. The uni- 
form decay of the scalar particle in the CM frame 
leads to a flat y e spectrum which is in contrast 
to the 1/i/g spectrum expected for the DIS NC 
background at fixed quark momentum fraction x. 
The DIS NC background for event topologies of 
type (B) with a "right" sign e is suppressed by 
exploiting both the angular distribution of the e 
and that of the highest P± jet found in the az- 
imuthal hemisphere opposite to the electron [jllj . 
Event topologies of type (B) with a "wrong" sign 
e are essentially background free. The misiden- 



Table 1 

Exclusion upper limits at 95% CL on A' lj7i , for Mq = 



150 GeV and M x o = 40 GeV for 7-like and Z-like X i- 
ha 

and are quoted at 95% CL. 



Indirect limits have been scaled to Mq = 150 GeV 



HI limits 



Indirect limits 



7-like Z-like Value Process [Rcf.] 



^'111 
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tified DIS CC or photoproduction background in 
event topologies of type (C) are suppressed by a 
combination of loose lepton (e or ^i's) identifica- 
tion and/or stringent f± requirements (f± > 15 
or 25 GeV). A good signal detection efficiency is 
obtained in each possible channel, e.g. ranging 
from 25% to 80% at M s = 150 GeV depending 
on the event topologies and model parameters. 

Except for a slight (~ 2.4er) excess around 
Mq = 70 GeV [|§jl| for "right" sign events 
in (B), an excellent agreement is found between 
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Figure 5. (top) Exclusion upper limits at 95% CL 
on X' 1 j 1 (regions above the curves excluded); (bottom) 
excluded domains in the M2 versus /x plane for a A' 
close to the upper limits at Mq = 150 GeV. 



observations and background expectations in all 
channels and exclusion limits at 95% CL are de- 
rived. The results for Ayj combining all con- 
tributing channels are shown in Fig. |5| (top) for 
M x o = 40 GeV. The existence of squarks with 
IjLp Yukawa couplings is excluded for masses 
up to - 220 GeV (240 GeV for M % o ~ 160 GeV) 
for coupling values X' 1 j 1 > \J^Ka em . The Fig. I 
(bottom) shows the excluded domains in the M2 
versus [i plane for Mq = 150 GeV and A' values 
near (or below) the sensitivity limit which is ob- 
tained at fixed M v o = 40 GeV. 
1 

To our knowledge, no direct searches in the 
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7? P -SUSY framework with A' 7^ have been per- 
formed yet by e + e~ and pp experiments. It is 
likely that in e + e~ collisions at LEP 200 where 
squarks are dominantly pair produced through 
gauge couplings, squarks masses up to near the 
kinematical limit of y/s e + e - /2 ~ 90 GeV can be 
reached within few years. In pp collisions at the 
TEVATRON, squarks can be produced in pairs 
or in association with gluinos. Given the cur- 
rent mass reach for scalar leptoquark and MSSM 
searches in DO and CDF experiments it is 
likely that Mg > 200 GeV can be probed through 
most of the MSSM parameter space even for 
Mg > Mq. Such a mass reach might be provided 
already by di-lepton data alone as was inferred in 
ref. g. 

From the analysis of the Ayj case involving 

the cLr and ul squarks, limits can be deduced 
on the Xijk by folding in the proper parton den- 
sities. Such limits are given in Table | at Mg = 
150 GeV. Also quoted in this table are the most 
severe indirect limits existing. The most stringent 
concerns A' ljfc with j = k and come either from 
the non-observation of neutrinoless double-beta 
decay (J = k = 1) or from constraints on the u e 
mass (j = k = 2,3). The limit from j3f3Qv decay 
depends on the gluino mass jl4j and is given here 
for Mg = 1 TeV. The limits on A' nfe from CC 
universality |15| affect only the d\. The limits 
inferred M] from M„ e constraints assume slep- 
ton and squark mass degeneracy and are valid 
only |l6| for A^ fc with j = k. 

In the presence of two simultaneously non- 
vanishing Yukawa couplings (e.g. A' lj7c and A^ fc 
with i =/= 1), resonant q production at HERA 
can be directly followed by a lepton flavor vio- 
lation (LFV) decay leading to fi + jet or r + jet 
signatures. Relevant analyses for these quasi- 
background free channels have been performed by 
the HI and ZEUS |o) collaborations. Exclu- 
sion limits in the context of ^ P -SUSY have been 
derived by ZEUS |o) based on L = 3.8 pb^ 1 of 
data and taking into account the branching ratio 
for gauge decays involving a pure 7 LSP (more 
unfavorable cases of mixed gaugino-higgsino mass 
eigenstates have not been considered) and assum- 
ing that only one type of squark is produced. The 
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Figure 6. Exclusion upper limits on X' 2 j k (top) 
and \' 3 j k (bottom) at 95% CL for squarks produced 
through \[ j k (regions above the curves excluded) . As- 
suming a pure 7 LSP, the d k limits are derived for 
^nfe — Kjk with i — 2,3 while the v? and t lim- 
its are derived for X'lji = ALj. with i = 2,3. The 
t limits assumes a mass mixing of cos 9 t = 0.5 and 
Mj < Mt- The upper (lower) solid curves are the 
results obtained for < M q (My ~ Mq). Indi- 
rect limits from other experiments |^lj are shown as 
dashed and dotted curves. 



results are shown in Fig. g. A sensitivity compa- 
rable or better than existing indirect LFV lim- 
its |2l| is obtained at Mq ~ 150 GeV for some 
coupling combinations. Squarks masses up to 
~ 210 GeV (depending on the LFV couplings) 
are excluded for A' > ^/Ana em . 

4. CONCLUSION 

Selectrons and squarks of the Minimal Su- 
persymmetric Standard Model were searched 
through pair production at HERA and masses up 
to 65 GeV were excluded in a new region at large 
M 2 and [i < 0. 

Squarks of i?-parity violating supersymmetry 
were searched through direct resonant production 
via Yukawa couplings \' ljk by the HI and ZEUS 
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experiments at HERA. 
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Assuming that one of the A^, dominates, three 
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families of event topologies were identified for R- 
parity violating and gauge decays of squarks. It 
was found that squark decays via a A' coupling 
into I + q final states dominate only at largest 
accessible masses, while elsewhere squarks un- 
dergo mainly gauge decays into a quark and a 
(generally) unstable gaugino-higgsino. No signif- 
icant evidence for the production of squarks was 
found in any of the channels and mass dependent 
limits on the couplings were derived. The exis- 
tence of first generation squarks with masses up 
to 220 -> 240 GeV (depending on the MSSM pa- 
rameter values) are excluded by HI at 95% con- 
fidence level for A'^ > ^47ra em . These limits 
extend beyond the current reach at other exist- 
ing colliders. At Mq = 150 GeV, the upper limits 
obtained for the four couplings X' ljk with j ^ 1 
and j ^ k are comparable or better than the most 
stringent indirect limits existing . 

Considering a combination \' x - k x A^ fc ^ with 
i ^ 1, squarks can undergo direct lepton flavor 
violation decays leading to /x + jet + X or r + 
jet + X signatures. For some coupling products, 
a sensitivity comparable or better than the most 
stringent indirect limit existing was obtained by 
ZEUS in the case of a "pure" 7 LSP. 
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